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abstract: Recent experiments show that an avalanche initiated from 
a point source propagates downwards by invading a triangular shaped re- 
gion. The opening angle of this triangle appears to reach 180° for a critical 
inclination of the pile, beyond which avalanches also propage upwards. We 
propose a simple interpretation of these observations, based on an extension 
of a phenomenological model for surface flows. 



Well-controlled avalanche experiments have recently been performed by 
Rajchenbach in (narrow) rotating drums [1], and by Daerr and Douady on 
a thin layer of sand atop a wide inclined plane (of high surface friction) at 
various angles of inclination [2]. In the latter experiments, an avalanche is 
triggered locally (e.g. by a probe) and its subsequent evolution is observed. 
One sees very clearly that below a certain inclination angle, the avalanche 
rapidly dies out, leaving only a small perturbed region around the initi- 
ation point. When the angle exceeds a first critical angle, the avalanche 
progresses downwards indefinitely; the perturbed region has a triangular 
shape with its apex at the initiation point. The opening angle of this tri- 
angle increases with the inclination and appears to reach 180 degrees at 
a second critical slope, beyond which the avalanche not only propagates 
downward, but also upward, progressively invading the whole sample. That 
avalanches can move upward has also been observed by Rajchenbach [1]. 
We want to show in this short note that these observations can be qualita- 
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tively (and perhaps quantitatively) understood within the so-called 'bcre' 
model that was introduced to describe surface flows and avalanches [3]. 

The starting point of the 'bcre' phenomenological description is to 
recognize that the simultaneous evolution of two physical quantities must 
be accounted for, namely: 

• the local 'height' of immobile particles, h(x, t) (which depends both 
on the horizontal coordinates x = (x,y) of the considered surface element 
and on time t) 

• the local density of rolling particles lZ(x,t), which can be thought of 
as the thickness of a flowing layer of grains. 

The presence of two variables, rather than one, crucially affects the 
"hydrodynamic" behaviour at large length- and time-scales and, as shown 
previously [3] can account for the Bagnold hysteresis effect within a simple 
one-dimensional treatment. (The Bagnold angle is the difference between 
the "angle of repose" of a surface after an avalanche, and the "maximum 
angle of stability" at which a new avalanche is initiated.) Variants of the 
model also allow segregation and stratification [4, 5] or the formation of 
sand ripples [6] to be understood in a simple way. 

In what follows, we extend the previous discussion of hysteresis [3] to 
2+1 dimensions (the case of a tilted planar surface) and show that the 
occurence of triangular avalanches can be understood. We also quantify 
the behaviour of backward propagating avalanche fronts and clarify the 
physical mechanisms involved in these. 

In constructing plausible equations governing the time evolution of the 
two quantities h and 1Z, we considered a regime where the rolling grains 
quickly reach a constant average velocity vq along the x direction (reflecting 
the balance between gravity and inelastic collisions with the immobile bed). 
We also assumed lZ(x, t) to be small enough in order to discard all effects 
of order 1Z 2 (for example, vq might depend on 1Z; we return to this point 
later). Thus we write [3]: 

an on n d 2 n n d 2 n ri _. 

-m = + ^ w +D ^W~ T [{nh {h}] (1) 

where the A)'s are (bare) diffusion constants, allowing for the velocity 
fluctuations of individual grains, and T describes the rate of conversion of 
rolling grains into immobile particles (or vice- versa). Correspondingly, the 
evolution of h(x, t) reads: 

^ = r [{K}, {h}} (2) 

since the total number of particles is conserved, and since by definition, 
the only mechanism by which the local number of immobile particles can 
change is by conversion into rolling particles. 
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Now, each rolling particle can, after colliding with the immobile bed, 
either come to rest or dislodge more particles. The rates at which these two 
processes occur obviously depend on the local geometry of the static grains 
near the surface; for simplicity we assume that this enters only through the 
local slope = —dh/dx, and the local curvature dO/dx. (By convention, 
we adopt 9 > for piles sloping downward in the positive x direction.) The 
probability of grains sticking is obviously a decreasing function of 9, while 
the probability for each grain to dislodge more wobbly particles increases 
with 6. Hence, for a certain critical value 9 = C (which we shall associate 
below with the angle of repose), the two effects on average compensate. For 
9 close to C , we thus expect: 

r [{ny, {h}\ = n U{e - e c ) + j'(e - e c ) 2 + ndo/dx..] (3) 

For 9 < C rolling grains, on average, disappear with time. On the contrary, 
for 9 > C , the rolling grain density proliferates exponentially, at least ini- 
tially (the nonlinear terms in TZ, neglected above, will then come into play). 
In the following, we shall in Eq.3 only retain the linear term in {0 — C ), 
although the quadratic term (7') can be important in some circumstances 
[7, 3]. The term in curvature, nd9/dx reflects the physical expectation that 
local 'humps' will tend to be eroded by a flux of rolling grains, while local 
'dips' tend to be filled in. Note that the coefficient 7 (which has dimen- 
sions of inverse time) can be interpreted as a characteristic frequency for 
collisions between rolling grains and the static substrate. 

So far, Eq. (3) assumes that the process by which rolling grains dislodge 
immobile grains is purely local. This might not be so - first of all, our con- 
tinuum description cannot be extended below the size of the grains (which 
we shall call a); the process by which a grain starting to roll destabilizes 
the grain which was just above it already leads to nonlocal terms in T. 
Furthermore, the momentum exchange between the rolling particle and the 
static bed can induce longer range effects through slight displacements of 
strings of contacts within the substrate. Mathematically, these non local ef- 
fects are described by adding higher order gradients (in TZ) to Y. However, 
the effect of the first two gradients in the expansion can be absorbed as 
a renormalisation of the values of vo, D || and D ± introduced in Eq. (1) 
above, to new values we shall denote v,D\\ and D±. 

Hence the value of D\\, which will turn out to play a crucial role in the 
following, reflects two separate effects: fluctuations of downhill velocity on 
the one hand, nonlocal dislodging effects on the other. The first contribution 
is of order Vq/j (recall that I/7 is a collision time), whereas the second is 
of order 7a 2 (or larger if the long range effects mentioned above turn out to 
be important). Both effects are thus of the same order of magnitude if, as is 
reasonable to expect, 7 ~ vo/a. We now ask what will happen to an initial 
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'pulse' of rolling grains. Will it progressively disappear with time, leaving 
the pile in a (globally) metastable state, or will it induce a 'catastrophic 
landslide'? As shown previously in 1+1 dimension [3], this depends on the 
initial angle of the slope 9{t = 0) = 9q. Suppose that the pulse was created 
at an arbitrary point which we choose as x = y = 0. After a time t, the 
density of rolling grains at site x is approximately given by [3]: 



lZ(x, t) = -j^=^^ exp 

4irJDnD ± t 



i(e -e c )t- ix ~ vt) 



AD\d ADj_t 



(4) 



where we neglect the modification of 9q brought about by the erosion pro- 
cess, and where the dependence on the transverse coordinate (y) has now 
been included. 

For #o < 9 C , the amplitude of the rolling grain pulse quicky reaches 
zero, after a time T ~ [t(©c ~~ ^o)] _1 - The length of the eroded region is 
thus finite, and diverges as i = v[y(@ c — ^o)] -1 . (Here and below, various 
logarithmic corrections to the leading algebraic behaviour are ignored.) If 
we now turn to slopes steeper than critical (9q > © c ) and look at the 
rolling grain density for a fixed X clS cl function of time, one sees that 
for 0o < @d = @c + ^ 2 /4-D||7, the rolling grain density grows, reaches a 
maximum, and then decreases with time. Assuming that x,y 3> a, the time 
t*(x,y) at which this maximum is reached is given by: 




and v = \Jv 2 — AD ^{9$ — 6 C ). We thus find there is a (spatially-varying) 
maximum value of the rolling grain density, given by lZ(x,t*). 

It is reasonable now to define the "avalanched" region, where grains 
have moved significantly, by the criterion 1Z(x,t*) > 7£ m i n , with TZ m i n a 
small threshold. The locus of points where the equality is satisfied defines 
the edges of the avalanche; to within logarithms it is easily shown to obey 



y 2 = (r + ^x^j - x 2 (6) 

where tq is a constant of order a (dependant on 7£ m i n ). This is a hyperbola, 
which asymptotes to two straight lines 
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as one moves away from the initiation point. Thus the model indeed predicts 
that the avalanching zone has (asymptotically) a triangular shape, with an 
opening angle if given by: 

K> - ^ ^ < 8 > 

Hence, if = when v = v , i.e. when 9q = C , and ip = tt when 6>o = &d- 
This analysis shows that tan(^) should diverge as (6^ — Oq)" 1 / 2 as the 
initial slope 6q approaches the second critical inclination angle, 0^. Quali- 
tatively at least, the model provides a simple explanation for the intriguing 
behaviour of triangular avalanches reported by Daerr and Douady [2]. 

We now remark on the propagation of backward avalanches for 6q > 0^, 
which can in fact be addressed within a purely one dimensional framework 
(dropping the y variable); our comments are relevant to Rajchenbach's 
experiments on narrow rotating drums [1]. If one looks at Eq. (4) at a 
given instant of time t, it tells us that the points where 1Z has reached a 
certain value 1Z* are given by: 

x±{t) = vt± 2^D||t[c + 7(0o - 6 c )t] ~ V±t (9) 

where the last inequality applies for t 2> 7 _1 and 

V± = v ± 2y/D in (0 -e c ) (10) 

The constant c$ depends on TZ*, but its precise value does not matter for 
large times (again we neglect logarithmic corrections). Interestingly, Eq. 
(10) means that for #o < @d, the two 'fronts' delimiting the zone where 
rolling grains are present at any given time are both propagating downhill 
(both have positive velocities). On the other hand, for #o > 0<i) one of 
the front moves with a velocity V- < 0. This is another way to say that a 
local perturbating pulse causes reorganization of the slope both uphill and 
downhill of where it started [3] ; the above argument quantifies the onset of 
this effect. The backward- moving front of dislodged grains is very clearly 
observed experimentally [1] and it would be interesting to test for the slope 
dependence of the propagation rate implied by Eq. (10). The experiments 
of Daerr and Douady [2] indeed suggest that the onset of the uphill-moving 
front might be directly identifiable as the maximum angle of stability, at 
least for thin inclined layers of sand; this interpretation of the 'spinodal 
angle' is exactly the one suggested in Ref.[3]. 

Note that one might have expected naively the upward-travelling front 
to propagate backwards only in a diffusive manner. The fact that it moves 
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with a finite velocity V_ is the result of the local diffusion constant D\\ al- 
lowing small backward motions, which are then amplified by the conversion 
effect 1 . It is important anyway to realize that, because the D\\ value reflects 
'nonlocal dislodgement' terms as well as the spread of downhill velocities 
(the latter is represented by the dispersion term A)||) the mechanism does 
not require that any individual grains are actually moving uphill. (On this 
point of interpretation we disagree with the remarks of de Gennes [8]. We 
also disagree with his view that the diffusive terms should anyway represent 
an ignorable correction to the leading hydrodynamic behaviour - for, even 
if very small, the diffusion constant controls a higher order derivative and 
can radically alter the long time behaviour [8].) 

Our model thus predicts the following scenario in the case where the 
slope is slightly steeper than (6q = ©^ + <5). When some external per- 
turbation creates a small local pulse of rolling grains, the 'rolling front' 
propagates downwards at velocity ~ v, and upwards with a smaller veloc- 
ity ~v5/ (@d — © c ) 2 - Once the downward front of grains hits the bottom of 
the silo, or of the rotating drum, the accumulation of grains creates there 
a bulge of immobile particles which, as pointed out in [3], moves upward 
at velocity 772. max . This occurs only if the time needed to create the bulge 
(L/v, where L is the linear size of the pile), is shorter than the time needed 
to complete the avalanche, which is given by L/772 max . Otherwise, the 
avalanche is already extinct when the bulge is created, so that backward 
propagation of the bulge (which requires nonzero TV) is precluded. A closely 
related mechanism is present for a silo being filled at a steady rate from a 
point source, and is involved in the stratification effect observed in grains 
of different sizes and roughness [5]. 

Note finally that within the model, for 9q > ©d 5 the rolling grain den- 
sity at any given point diverges to infinity at long times. Obviously, this 
divergence is unrealistic; but it can easily be corrected if higher powers of 7Z 
are included in Eq. (1). The simplest assumption 3 is to write 7 = 70 — 71 72. 
, representing the fact that for larger values of 1Z, not all the rolling grains 
interact with the solid phase. (This leads to an effective reduction of 7.) 
This correction term will cause the growth of 1Z to saturate at some limiting 
value 72 max = 70/71; the resulting non linear equation is known as the kpp 
equation (Kolmogorov, Petrovsky and Piscounov), and sustains propagat- 
ing fronts between the 'flowing' region and the region at rest. The velocity 

1 The same remark applies to the transverse propagation effect, leading to triangular 
avalanches, which is mediated by D± 

2 Note that since O d — O c is of the order of a few degrees, this upward velocity is 
actually expected to be of the same order of magnitude as v, unless S can be made much 
smaller than 1° 

3 De Gennes has alternatively suggested to describe the saturation effect by substitut- 
ing -yTZa/(JZ + a) to -yTZ [9]. 
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of these fronts depends on their shape, but its order of magnitude is still 
given by the above values of V± [10]. The main features of the hysteresis 
scenario, including triangular and backward avalanches, as described above 
is thus qualitatively unaffected. 
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